Tunneling Time of Bose-Einstein Condensates on Real Time Stochastic Approach 
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We study tunneling processes of Bose-Einstein condensate (BEC) on the real time stochastic 
approach and reveal some properties of their tunneling time. An important result is that the 
tunneling time decreases as the repulsive interatomic interaction becomes stronger. Furthermore, 
the tunneling time in a strong interaction region is not much affected by the potential height and 
is represented by an almost constant function. We also obtain the other related times such as the 
hesitating and interaction ones and investigate their dependence on the interaction strength. Finally, 
we calculate the mean arrival time of BEC wave packet and show the large displacement of its peak 
position. 

PACS numbers: 03.75.Lm,42.25.Bs,03.65.Xp,03.65.Sq 



CZ2 

wo; 

•4-* 



-a 
c 

o 
o 



> 

(N 
O 
oo 
O 

o 



X 



I. INTRODUCTION 

The Bose-Einstein condensates (BECs) of trapped 
atoms, first realized in 1995 [H-Q, are ideal systems for 
studying macroscopic quantum phenomena. This is be- 
cause the systems are dilute, weakly interacting ones and 
we can easily control the configuration of the trap and 
even the strength of interatomic interaction. The BEC 
is described by a macroscopic wavefunction with nonlin- 
ear effects due to interatomic interactions and show the 
macroscopic tunneling phenomena such as Josephson os- 
cillations. Many studies on BEC tunneling through vari- 
ous kinds of potentials |4H6j have already been done. We 
expect that the BECs give us new insights for nonlinear 
dynamics and macroscopic tunneling phenomena. 

The problem of tunneling time for quantum particle is 
controversial, mainly because the time is not represented 
by an operator but is a parameter in quantum mechan- 
ics. It has been investigated for decades and many defi- 
nitions of the tunneling time have been proposed |7Hl8j. 
For example, the phase time Q is expressed by an energy 
derivative of phase shift of transmitted wave from a po- 
tential barrier. The Larmor time (8rll0| is obtained from 
the Larmor precession angle caused by a magnetic field in 
the potential barrier region. The Biittiker-Landauer time 
|1TH13| is defined from transmission coefficient through 
a barrier with a small oscillation of height. There are 
methods based on the concept of "particle path" such 
as Bohmian mechanics [3, [II], Feynman path integral 
[Trj , Nelson's stochastic approach [l8rl2l] | and so on. 
The cold atomic gas systems are promising to investigate 
the tunneling time of quantum particle, since their dy- 
namical behaviors can be observed in detail and one can 
even control several parameters such as the configuration 
of trap, the strength of interatomic interaction and the 
internal degree of freedom such as spin. 
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In this paper, we focus on the tunneling time asso- 
ciated with BEC wave function and its dependence on 
the interaction strength. But most of the existing def- 
initions of the tunneling time is basically based on the 
Schrodinger equation for one-particle tunneling phenom- 
ena. The BEC dynamics is well described, particularly 
near the zero temperature, by the Gross-Pitaevskii (GP) 
equation, which is the non-linear Schrodinger equation. 
We can not adopt the previous definitions of the tun- 
neling time to the BEC wave in a straightforward way. 
Zhang and his co-worker estimated the tunneling time of 
the BEC wave packet qualitatively form its peak motion 
[22j ■ Their results show that the tunneling time strongly 
depends on the strength of the non-linear interaction. 
But the appearance of the negative tunneling time is in- 
evitable in their approach. 

For our purpose, we use the Nelson's stochastic ap- 
proach. In Nelson's stochastic quantum mechanics [2a , 
|24| . the random variable x(t) representing the motion 
of a quantum particle is subject to the stochastic dif- 
ferential equations. It gives the real time trajectories of 
the quantum particle as sample paths. From ensemble of 
the sample paths, we can reproduce the predictions given 
by the ordinary quantum mechanics. Since the Nelson's 
stochastic approach utilizes the real time trajectories, we 
acquire direct information about the tunneling time. We 
note that it can be extended to the system described by 
the non-linear Schrodinger equation [25rl27j . 

In this paper, we explicitly calculate the times related 
to tunneling phenomena of the BEC wave packets by 
means of the Nelson's approach. Their dependence on 
the interaction strength is an interesting subject. We 
note that the GP equation can also be derived by ap- 
plying the mean field approach directly within the for- 
mulation of the Nelson's quantum mechanics, as will be 
shown in this paper. We perform numerical calculations 
for a quasi one-dimensional system with a rectangular po- 
tential barrier, that is, one with a barrier in x-direction 
and a strong confining harmonic potential in yz-plane. It 
will be seen from the accumulated sample paths how the 
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tunneling time of BEC depends on the non-linear inter- 
action strength: it decreases as the interaction becomes 
stronger. Furthermore, it becomes constant in a strong 
interaction region. To explain those results, we analyze 
a simple model with a double well potential. Next, we 
focus on the hesitating behaviors found in the sample 
paths, which are due to a strong interference between the 
incident and reflecting waves. We evaluate the hesitating 
time and find that it is much affected by the interaction 
strength. From the calculations of the tunneling and hes- 
itating times, we confirm that the Hartman effect Q is 
violated when the non-linear interaction is switched on. 
Finally, we also calculate the arrival time distribution of 
the BEC wave packet and the mean arrival time from it. 
The result suggests a large acceleration in the motion of 
the BEC wave packet in the presence of a potential bar- 
rier and the non-linear interaction. Nelson's approach 
can reproduce the arrival time distribution. We argue 
that the displacement of the peak position accounts for 
the "acceleration". 

This paper is organized as follows. In Sec. II, we re- 
view the Nelson's quantum mechanics and derive the CP 
equation directly in the stochastic approach, using the 
mean field approach. In Sec. Ill, we numerically calcu- 
late the tunneling time in the stochastic approach. In 
Sec. IV, we discuss the mean arrival time of the BEC 
wave packet. Section V is devoted to summary. 



II. REAL TIME STOCHASTIC PROCESS AND 
GROSS-PITAEVSKII EQUATION 

In this section, we briefly review the Nelson's stochastic 
quantum mechanics with N degrees of freedom [23l . [2~3 . 
1271 ] and show its equivalence to the Schrodinger equation. 
Next, applying the mean field approach to the dynamical 
and kinematical equations, we derive the CP equation. 

The first assumption of the Nelson's stochastic quan- 
tization is to set up the stochastic differential equations 
for the particle position Xj(t) for i = 1, • • ■ , N as 

dxJt) = b(X(t),t)dt + J— dWi, (1) 
V rn t 

dxi{t) — Xi(t + dt) — Xi(t) , 
for forward time evolution and 

dxAt) = bJX(t),t)dt + J— dWi*. (2) 

V 

dxi(t) — Xi(t) — Xi(t — dt) , 

for backward one with dt > 0, particle masses rrii and 
the notation, X(t) = (xi(t), ■ ■ ■ ,X]y(t)) . Here Wi and 
Wi* are three-dimensional independent standard Wiener 
processes 

E[dW Ll (t)} = E[dW L „(t)} = , (3) 
E[dW hti (t)dW l2j (t)} = E[dW* lui (t)dW* l2 j(t)} 
= 5i 1 i 2 5 ij dt, (4) 



where !?[•••] means an ensemble average and indices I 
and i represent I — x,y,z and i = 1, • • • , N , respectively. 

The second assumption is the Newton-Nelson's equa- 
tion of motion as follows: 



i a i {t) = -V i V{X(t),t) 



(5) 



where V(X(t),t) represents the sum of the external po- 
tential V ex (xi,t) and interaction potential Vi nt (xi — Xj) 
as V = Ei V m (x it t) + (1/2) YUftj Vtatfa - x,) . The 
mean acceleration a,i(t) is defined as 



«*(*) = \{DD m +D m D)xi(t), 
with the mean forward time derivative D , 
-f(X(t + dt))-f(X(t)) 



Df(t) = lim E 



dt 



and the mean backward one , 

-f(X(t))-f(X(t-dt)) 



D*f(i) = lim E 



dt 



X{t) 



X(t) 



(6) 



(7) 



,(8) 



where E[- ■ ■ \X(t)} means the conditional expectation. 
Let us define the current and osmotic velocities as v = 
(b + 6*)/2 and u = (b - 6*)/2. Then, the Newton- 
Nelson's equation of motion ([5]) becomes the equation 
for the current velocity: 



d 



E 



■2 

2nij 1 



1 



(9) 



The stochastic processes of the random variables Xi (t) 
in Eqs.([T]) and @ are equivalently formulated by means 
of the distribution function P(X,t) which satisfies the 
Fokker-Planck equations, 



dt 

i 

for forward time and 
d 



dt 



(10) 



(ii) 



for backward one. It is known that the forward and 
backward velocities have a certain relation bi — b*i = 
^VilnP [24]. From Eqs.lJISJ), (HI]) and m = = 

J^— Vj In P , one derives the equation for the osmotic ve- 
locity as 

d_ 

ot 



Ui 



E 



2m j 



(12) 



3 



Equations ([9]) and (fT2j) are called the dynamical and kine- 
matical equations, respectively, and their combination 
can be transformed into the Schrodinger equation, 



£^Vf + U(X,t) W(X,i)(13) 
L — ' 2m, ' 



with the substitutions, Vi = -^Im[V, In and Ui = 
— Re[Vj In 1 !'] . The probability of the particle posi- 
tion P(X,t) corresponds to the absolute square of the 
Schrodinger amplitude P(X,t) = \^{X,t)\ 2 . Thus, the 
stochastic differential equations ([T]) and ((2]) are equiva- 
lent to the Schrodinger equation. 

Next, we apply the mean field approach to the Nel- 
son's stochastic quantum mechanics and derive the GP 
equation. We average the i-th dynamical and kinematical 
equations over the variables (x±, • • • , Xi—i,Xi+i, ■ ■ ■ ,xn) 
as 



dt 



= 5>« 



— E x 



~ — Ef 

dt 



2m 
1 



-ViV 



in 



Ui 



2m 



(14) 



(15) 



where E Xi [■ ■ ■ ] means the conditional expectation as 

E X .[F(X)] = jY[dx 3 F(X) Pl (X,t), (16) 



P(X,t) 



m^d Xj p(x,t) 



(17) 



Here, we apply the mean field ansatz and put the fac- 
torized probability distribution for the particles position 
as P(X,t) = Ylj Pj(xj,t) . Since a BEC phase is under 
consideration, every particle should have the same prob- 
ability distribution P\{x,t) = ■ ■ ■ = Pjy(x,t) — P(x,t) . 
Then, all the current and osmotic velocities should be 
in the same forms irrespective of the index i, Vi(X,t) — 
v(xi,t) = Vi and Ui(X, t) = u(xi, t) = Ui, and Eqs. (fT4")) 
and ([TBI become 
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Vi(ui ■ Vi) 



(18) 
(19) 



Now, we take the contact-type interaction potential, 

V in t(xi - Xj) = gS(xi - Xj) (20) 



where g = 4irh 2 a s /m with s-wave scattering length 
a s . Note that a repulsive interaction g > is con- 
sidered for BEC. Utilizing the transformations €>j = 
^Im[Vi ln^>] and iii = ^Re[V.; Inip] , we obtain the non- 
linear Schrodinger equation from Eqs. (IT8l) and (|T9)) as 



.8.. . 
th—ip(x,t) = 



2m 



■g(N-l)W 



(21) 



With TV ~ N - 1 and * = VNip, it is the GP equation. 
The corresponding stochastic differential equations are 
given by 



dx(t) = [v(x(t),t) + u(x(t), t)} dt + \ —dW , (22) 



dx(t) = \v(x(t),t) - u(x(t), t)} dt + \ —dW* (23) 

" m 



III. TUNNELING TIME OF BOSE-EINSTEIN 
CONDENSATE 

We consider a system with an external potential in 
x-direction, V ex (x) , and a strong confining harmonic po- 
tential in j/z-plane. Suppose that the transverse confine- 
ment is very strong and the transverse excitations are 
forbidden. Then the total wave function ty(x,t) can be 
approximated as 



V(x,t) ~ i/)(x,t)<l>(r ± )t 



-iojj_t 



(24) 



where %p(x,t) is the wave function in x-direction mode 
and <f>(r±) (rj_ — \J y 2 + z 2 ) is that of the trans- 
verse ground state. This system is called a quasi one- 
dimensional system, and we perform numerical calcula- 
tions with a rectangular barrier potential for V cx (x) . 

For the weak interaction, (j){ r ±) can be replaced by the 
Gaussian form 



-) = 



■ exp 



mu>±rj 
2~H~ 



(25) 



Then, the GP equation for ip(x, t) is given by 
d 2 



dx 2 



+ V ex {x)+g\^{x,t)Y 



,(26) 



with a; in a unit of ai = y/h/muj±, t in a unit of wj 1 , the 
dimensionless wave function ijj/a±_ — ¥ ip , and the nonlin- 
ear interaction 2a s N/a± — > g . We consider a rectangular 
barrier potential as 

, x < (region I) , 

V cx (x) = { V , < x < d (region II) , (27) 
, d < x (region III) , 
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with the unit of hujj_ . We assume that the initial wave 
packet with a variance (Ax) 2 has the Gaussian form as 



i 



2?r(Ax) 



1/4 



exp 



(a: - x ) s 
' 4(Ax) 2 



+ik (x — Xq) 



(28) 



where xq and fco are the initial center position and mo- 
mentum of the wave packet, respectively. 
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FIG. 1: (Color Online) (a) The transmission rate of the 
BEC wave packet versus the potential width for the inter- 
action strength g — 0, 1, 3, 5 . (b) The center momentum of 
the transmitted wave packet versus the potential width for 
g = 0,1,3, 5. The other parameters are Vb = 1.4(fcg/2) and 
k = 0.6. 



energy into the kinetic one particularly in high-density 
regions. The energy conversion is also responsible for the 
dependence of the center momentum 



fc = dkk\ip T (k)\' 



(30) 



on g, as is shwon in FigQ]-(b) . Finally, we refer to the 
initial position dependence of the transmission rate. The 
interaction energy is also being transformed into the ki- 
netic energy even outside of the potential barrier and 
broadens the width of the wave packet. It affects the 
transmission rate in Fig J5] for xo = —100,-150,-200, 
but the dependence is not significant in our choice of the 
parameters. So we show only the results for the initial 
position xq = —150 below. 
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FIG. 3: Typical transmission sample path with hesitation. 
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FIG. 2: (Color Online) The transmission rate of the BEC 
wave packet versus the potential width for the initial position 
xo = —100,-150,-200. The other parameters are g = 3 , 
Vo = 1.4(fcg/2) and fc =0.6. 

We now investigate the tunneling of a BEC wave 
packet. In Fig[T]-(a), we calculate its transmission rate 
as 

/>oo 

T= / dx|V> T 0Mf)| 2 , (29) 
Jd 

where ipr is the transmitted wave packet and t{ is the fi- 
nal time of the scattering problem. We can find that the 
transmission rates depend very much on the nonlinear in- 
teraction strength g . This dependence is mainly because 
the repulsive interaction converts part of the interaction 



Let us turn to calculate the tunneling time of the BEC 
wave packet. Figure [3] shows a typical transmission sam- 
ple path, calculated by Eg. (|2"31 with "backward time evo- 
lution method" 0, 0j| • The tunneling time tt is defined 
as the averaged time interval in which the random vari- 
ables x(t) stay in the barrier region II, 

^Sample t/ 

tt = Yl m / QMW' (3i) 

; JV Samplo Jt t =0 

( 0, x < 

8(x) = < 1 , < x < d (32) 
[ 0, d < x, 

where i$ = and tf are the initial and final times for 
the scattering problem, respectively. As shown in Fig|3l 
the random variable x(t) of the transmission sample path 
stays in front of the potential barrier [l8j |. This hesitat- 
ing phenomenon is due to a strong interference between 
the incident and reflecting waves. We define the hesitat- 
ing time th as the averaged time interval in which the 
random variables x(t) pass through some region in front 
of the potential barrier (— dn < x < 0) . The interaction 
time T\ is also defined as a sum of the tunneling time 
tt and the hesitating one rjj . Thus the interaction time 
tj = tt + tjj represents the passage time through the 
region (— dn < x < d) . 
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solutions are given by 

9 =o — 

2=3-- ; AN = -AN cos(ut) , (35) 

- AiVo-^sin(^) , (36) 

with the frequency hw = y/2J(NU + 2 J) . Since we are 
j | not interested in an oscillation of the two BECs here but 
only in the tunneling from the left to the right, we con- 
sider only for uit <§; 1 and ignore the order 0(t 2 ). Then, 
we obtain 



FIG. 4: (Color Online) (a) The tunneling time of the BEC 
wave packet versus the interaction strength g with the param- 
eters d = 7, Vo = 1.4(Jfeg/2) and k = 0.6. (b) The tunnel- 
ing time of the BEC wave packet versus the potential height 
Vo/(feo/2) for g = 0, 1,3, 5. The other parameters are d — 7 
and fco = 0.6 . 



We study how the interaction strength g affects the 
tunneling time. Figure 0J- (a) shows a the behavior of the 
tunneling time versus g. We can see that the tunneling 
time decreases as g increases. The nonlinear repulsive in- 
teraction accelerates the motion of the BEC wave packet 
in the potential region. Figure HJ-(b) shows the behavior 
of the tunneling time versus the potential height. The 
tunneling time with 5 = increases first as the potential 
becomes higher, which can be understood intuitively, but 
shifts to a decrease for the high potential. The latter be- 
havior can be explained by the Biittiker-Landauer time 
t = radj (h^2(Vo — E)) , which is also obtained by the 
Nelson stochastic approach for high and wide potential 
barrier with g = (18| . On the other hand, the tunneling 
time does not vary much with the potential height Vq in 
case of the strong interaction, and approaches to a con- 
stant value. Our results imply that the tunneling time of 
the BEC wave packet mainly depends on the interaction 
strength g, but not on the potential height. 

To explain the Vo independence of the tunnneling time 
for large g, we take a simple model of BECs in a double 
well potential. We suppose that two BECs are initially in 
a stationary state with the condensate particle number 
difference AiV(0) = and the phase difference Acf> = . 
Then, we add AN particles to the left well, AN(0) = 
—AN , and this number difference induces the tunnel 
current form the left well to the right one. Then the 
motion of the BECs between the two wells are described 
by the simultaneous equations [28] 



j t AN(t) = 2^A<p(t), (33) 
±A<f>(t) = -(u+?p)AN(t), (34) 



where N, J, and U represent the total condensate parti- 
cle number, the tunnel coefficient, and the interaction 
constant, respectively. It is assumed in derivation of 
Eqs.flSS]) and (jUJ) that AN /N and A0 are small. Their 



Ac/) 



AN hu 2 t 
~lv 2J~ 



' ANq 2Jt 
N h 

ANq NUt 
. N h 



(for U < J/N) , 
(for U > J/N) . 



(37) 



The phase difference A<j) corresponds to the tunneling 
current velocity, so larger A(j) implies smaller tunneling 
time. One notices the following two features of the phase 
difference A(f>. First, it increases monotonically as U 
does. Second, while it depends on J in the weak in- 
teraction limit U <C J/N, it becomes independent in the 
strong interaction limit U ^ J/N. As U and J in the 
toy model with the double well can be identified with g 
and the potential height Vo in the model of the potential 
barrier, respectively, the arguments just above explain 
the behaviors of the tunneling time in FigJ3]- 




FIG. 5: (Color Online) (a) The tunneling time of the BEC 
wave packet versus the potential width d for g — 0,1,3,5. 
(b) The hesitating time of the BEC wave packet versus d 
for dn — 10 and g = 0,1,3,5. The other parameters are 
Vo = 1.4(fco/2) and k = 0.6 . 
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FIG. 6: (Color Online) The interaction time of the BEC wave 
packet versus the potential width d for g — 0,1,3,5. The 
other parameters are Vo = 1.4(/cq/2) and fco = 0.6 . 



Finally, we investigate the properties of the hesitat- 
ing and interaction times versus the potential width d , 
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which are shown in Fig|5j For g = 0, while the tun- 
neling time increases monotonically, the hesitating time 
decreases for large d. This is because that the particle 
of the transmission sample path for thick potential tends 
to have higher velocity and therefore to pass through the 
region in front of the potential barrier faster [29]. As g 
becomes larger, the tunneling and hesitating times de- 
crease. One can see that the hesitating time is affected 
by the interaction strength much more than the tunnel- 
ing time, as in Fig 03 This result can be explained as 
follows: The incident and reflecting waves make a strong 
interference in front of the potential barrier and create 
the high density region. There the non-linear repulsive 
interaction term g\ip\ 2 contributes strongly to the behav- 
ior of the sample path. Next, we refer to the dependence 
of the interaction time on d, as in Fig|5]. It is predicted, 
based on the method of phase time Q , that the "tunnel- 
ing time" for thick-enough barrier becomes independent 
of the barrier length for non-interacting systems, which 
is known as the Hartman effect. In the Nelson's stochas- 
tic approach, the tunneling time with g = grows but 
the hesitating time decreases, and their sum, the inter- 
action time, seems to approach to a constant value, as 
in Fig|6] This corresponds to the Hartman effect in Nel- 
son's stochastic approach. We remark that the Hartman 
effect is apparently violated in the presence of the non- 
linear interaction, g ^ . The violation of the Hartman 
effect for non-linear interaction has also been pointed out 
in Ref. (H . 



IV. ARRIVAL TIME OF THE BEC WAVE 
PACKET 

As seen in FigQ~|-(b) , the center momentum for the 
transmitted wave packet becomes larger than that for the 
incident one. In tunneling process, the non-linear inter- 
action reduces the tunneling and hesitating times. The 
acceleration of the quantum particle in the presence of 
the potential barrier and nonlinear interaction accounts 
for these results. In this section, we investigate the ac- 
celeration of the wave packet in view of the mean arrival 
time [11]. 

Introduce the arrival time distribution at the position 

x as 



P x (t) 



f*?dtMx,t)\i 



(38) 



and the mean arrival time T x as 

T x = [' dttP x (t) 



(39) 



One can calculate the difference between the mean arrival 
times with and without a potential barrier, 



Due to the nonlinear interaction, the center of the mo- 
mentum for the transmitted wave packet becomes large 
in Fig[T|-(b) and the mean arrival time T x should reflect 
this effect. In order to study the acceleration of the BEC 
wave packet in the potential barrier, we focus on the 
mean arrival time at the potential barrier edge x = d . 




FIG. 7: (Color Online) (a) The difference between the mean 
arrival times of the BEC wave packets with and without a 
potential barrier versus the potential width d for g = 0, 1, 3, 5 . 
(b) The difference between the interaction times of the BEC 
wave packets with and without a potential barrier versus the 
potential width d for g = 0, 1, 3, 5 . The other parameters are 
V = 1.4(fco/2) and k = 0.6. 

The results of AT x= d are shown in Fig|7]-(a) . At first, 
we see that AT X=C [ can become negative and that the 
non-linear interaction gives rise to a large | A.T , a;=t i | , which 
strongly suggests a big acceleration inside the potential 
barrier region. For comparison, we also show the dif- 
ference between the interaction times with a potential 
barrier tt and without it rf 00 in FigEKb) • The differ- 
ence Ati = t\ — Tj e becomes small as the interaction 
strength goes up, but does not become negative in con- 
trast to AT x= d- It indicates that the velocity of the trans- 
mitted wave packet does not exceed that of the free wave 
packet. Although the above results sound paradoxical at 
first glance, the Nelson's method can reproduce physical 
quantities in quantum mechanics and the arrival time dis- 
tribution can actually be obtained from the transmitted 
sample paths (FigJS]). 
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FIG. 8: (Color Online) (a) The arrival time distribution of 
the BEC wave packet P x= io(t). (b) The histogram obtained 
from an ensemble of transmitted sample paths at x — 10 . 
The other parameters are g = 5, d = 10, Vo = 1.4(fcg/2) and 
k = 0.6. 
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FIG. 9: (Color Online) (a) The absolute square of the am- 
plitude ip(x,t) for t = 120,160,200. (b) The probability 
distribution obtained from ensembles of the transmitted and 
reflected sample paths for t = 120, 160, 200 . (c) The prob- 
ability distributions of the transmitted and reflected compo- 
nents, which are obtained from ensembles of the transmitted 
and reflected sample paths for t = 120, 160, 200 . The other 
parameters are g = 5, d = 10, Vb = 1.4(fcjj/2) and ko = 0.6 . 

The key to understand these results is a displacement 
of the peak position of the wave packet. Figure 0(a) 
shows the wave packets at t = 120, 160,200. It tells us 
that before the peak of the incident wave packets reaches 
the potential barrier x = the peak of the transmitted 
one appears at t = 200 . Furthermore, only the front 
part of the incident wave packet seems to contribute to 
the transmission. This situation becomes clear in the 
Nelson's stochastic interpretation. From the transmitted 
and reflected sample paths, we can construct the prob- 
ability distributions for the respective components. Fig- 
ure [9] shows that the transmitted wave packet is con- 
structed mainly from the front part of the incident wave 
packet. As a result, a displacement of the peak posi- 
tion of the wave packet occurs and it looks like its large 
acceleration. This mechanism is similar to the superlu- 
minal tunneling of the photon [3(3, HH . We point out the 
strong dependence of the displacement on the interaction 
strength. 



V. SUMMARY 

In this paper, we have investigated the times related to 
the tunnneling of the BEC wave packet in the Nelson's 
stochastic approach. There are the three times, namely 



the tunneling time for which a particle is in a potential 
barrier, the hesitating one for which it stays in front of 
the barrier and the interaction one, give by a sum of 
the tunnneling and hesitating times. Applying the mean 
field approach to the Nelson's stochastic formulation, we 
derive the GP equation directly. 

According to the stochastic formulation, we have per- 
formed numerical calculations. First, it is found that the 
tunneling time decreases as the interatomic repulsive in- 
teraction becomes stronger and is not affected so much 
by the potential barrier height for the strong interaction. 
The dependence of the hesitating and interaction times 
on the parameters, especially on the interaction strength, 
has been revealed. It is seen that the Hartman effect is 
violated when the non-linear interaction is switched on. 

We have also calculated the mean arrival time of the 
BEC wave packet and have seen a large displacement 
of its peak position. This result implies that it is not 
adequate to define the tunneling time by the peak motion 
( or the group velocity ) of the BEC wave packet. 
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